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Abstract 

In the present note we study absolute linear Harbourne constants. These 
are invariants which were introduced in [2j in order to relate the lower bounds 
on the selfintersection of negative curves on birationally equivalent surfaces to 
the complexity of the birational map between them. We provide various lower 
and upper bounds on Harbourne constants and give their values for the number 
of lines s of the form p 2r +p r +1 for any prime number p and also for all values 
of s up to 31. This extends considerably results of the third author obtained 
earlier in [El- 
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1 Introduction 

Arrangements of lines were introduced to algebraic geometry by Hirzebruch in his 
papers concerning the geography of surfaces (i.e. construction of surfaces X with 
prefixed invariants cy(X) and 02 (A)), see [7], (T . 

Multiplier ideals defined by arrangements of lines were studied by Teitler El 
and Mustafa [8j. 

Recently arrangements of lines appeared in the ideas revolving around the Boun¬ 
ded Negativity Conjecture (BNC for short), see [3J for the background of the Con¬ 
jecture and 0 , m for the role of configurations of lines. Whereas BNC is relevant 
only over a field of characteristic zero, some related problems are of interest over 
arbitrary fields. In [2J the authors introduced and began to study linear Harbourne 
constants. These are certain invariants computed by configurations of lines in the 
projective plane. Even though the Bounded Negativity fails in positive characteris¬ 
tic, it is clear from Definition o that for a fixed d , the linear Harbourne constant 
H(d) is a finite number (because the number of combinatorial possibilities for in¬ 
variants of a configuration of d lines is finite). It is interesting to estimate these 
numbers because in particular they measure the discrepancy between combinatorial 
data sets, see [3] and those sets which come from geometric configurations defined 
over some fields. 

For the purpose of this note, a configuration L is a finite set of mutually distinct 
lines £ = {L \,..., L ^}. Given a configuration £, we define its singular set CP(£) = 
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{Pi,..., P s } as a set of points where two or more lines intersect. This is the same 
as the singular locus of the divisor Li + ... + L For a point P € T(£), we denote 
by mc(P) its multiplicity, i.e. the number of lines which pass through P. We have 
the following definitions. 


Definition 1.1. The linear Harbourne constant of a configuration of lines L in the 
projective plane P 2 (IC) is the rational number 


H(K,L) = 


d 2 -Efc=i m -c( p fc ) 2 


(1) 


The linear Harbourne constant of d lines over K is defined as the minimum 


H{K,d) := min H(K,£>) 


taken over all configurations £ of d lines. 

Finally the absolute linear Harbourne constant of d lines is the minimum 


H(d) := min H (K, d) 


taken over all fields 1C. 

In order to alleviate the notation we define first the set 


Q = {q = p r , p is prime, r € Z>o}. 

For an integer d, we define q(d ) as the least number q € Q satisfying 

d ^ q 2 + q + 1 

and r(d) as the largest number r € Q satisfying 

r 2 + r + 1 ^ d. 


Systematic investigation of absolute linear Harbourne constants H(d) was initiated 
in m- Results stated there and computer supported experiments have led us to 
formulate the following conjecture. 

Conjecture 1.2. For d ^ 2 let q = q(d') and let i := q 2 + q + 1 — d. 

If i ^ 2q — 2, then 

H(d) = h(d) 

where 

_ q 2 + q + l- i- £\(i)mi{i) - e 2 (i)m. 2 (i) - t q _i{i)(g - 1) - t q (i)g - t g+ i(i)(q + 1) 

£l(i) +£2(2) +tq- 1 (i) +t q (i) +t q+ i(i) 


with 


ei(i) 


m i(0 = q + 1 — i, m 2 (i) = 2q + 1 — i 


1 for 0 ^ i ^ q — 1 
0 otherwise 


( 1 for i > q + 1 
{ 0 otherwise ’ 


qi — q 2 — q 
0 


*5-1 (*) 


for i > q + 1 
otherwise 
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*«(») = 


qi 


2 q l - (* - 2)g - 1 


*5+1 (*) = 


q z + q — iq 
0 


Moreover for i = 2q — 1 we have 


for i ^ g + 1 
/or i > q + 1 

/or i ^ q + 1 
otherwise 


H(d ) = - 


<r - q 2 + 2g - 2 
< 7 2 + q - 1 


Remark 1.3. We do not know what happens for d such that d ^ (q(d) — l) 2 + q(d). 
The first such d is d = 32 with <7(32) = 7. See the end of the last section. 

This conjecture has been verified for d ^ 10 in [10| . In the present paper, we 
extend the range of the validity of the Conjecture to d ^ 31. This is our first main 
result. We repeat the results from [10] for completeness. 

Theorem 1.4 (Values of absolute linear Harbourne constants). For 2 ^ d ^ 31 we 
have 


d 

H{d) 



2 

0 


d 

H{d) 

3 

-1 

d 

H(d) 

22 

-108/29 « -3.724 

4 

-4/3 « -1.333 

14 

-54/19 « -2.842 

23 

-115/30 » -3.833 

5 

-3/2 = -1.5 

15 

-3 

24 

-4 

6 

-12/7 w -1.714 

16 

-16/5 = -3.2 

25 

-125/30 « -4.166 

7 

-2 

17 

-67/20 = -3.35 

26 

-129/30 = 4.3 

8 

-2 

18 

-24/7 ss -3.428 

27 

-135/31 fa -4.354 

9 

-9/4 = -2.25 

19 

-76/21 » -3.619 

28 

-140/31 ta -4.516 

10 

-29/12 » -2.416 

20 

-80/21 « -3.809 

29 

-145/31 ta -4.677 

11 

-33/13 « -2.538 

21 

-4 

30 

-150/31 ta -4.838 

12 

-36/13 « -2.769 


31 

-5 

13 

-3 




Table 1 Values of H(d) for up to 31 lines 


Remark 1.5. It comes as a surprise that the function H(d) is not decreasing with 
d increasing. 

We prove the Conjecture for any d = q(d) 2 + q(d) + 1, see Corollary [U 
The last assertion is a consequence of the following more general result. 

Theorem 1.6 (Lower bound on Harbourne constants). For d ^ 6 we have 

H(d) Js ——V4d — 3 + —. 

For d = q 2 + q + 1 with q € Q we have the equality. In this case H(d) = —q is 
computed by the configuration consisting of all lines in the finite projective plane 
P 2 (F ? ). 
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Theorem 1.7 (Upper bound for Harbourne constants). For d ^ 7 and mf/i r = 
r(d), we have 

r 4 + r 3_ y ,_ ( d _ 1)2 
r 4 + 2 r 3 -r-d 2 + d- 2 ' 

We will considerably improve this bound for some d in Proposition l5.il In order 
to prove Theorem 11.41 we introduce some new tools, which might be of independent 
interest in other areas of combinatorics and geometry. We discuss also how our prob¬ 
lem is related to the classical geometric problem on the existence of projective planes 
with certain numbers of points. Our investigations are accompanied by Singular [ 6 ] 
computations. The complete script of our program is provided in the Appendix. 


2 Initial data 

Let L = {L \,..., Ld) be a configuration of lines in the projective plane P 2 (IC). Let p. 
be the number of points where exactly k lines intersect. Then we have the following 
basic combinatorial equality 



( 2 ) 


Note that using this notation and taking into account ([5]) we can simplify the way 
H(K,L) is expressed: 


H , K £) = - Efc= 2 = d - Efc =1 m L {P k ) 

s s 


(3) 


Now our approach to computing or bounding Harbourne constants is based on 
the following idea. For a fixed d we consider the set 7 of all integral solutions 
T = (^ 2 , £ 3 , • • •, td) of the equality © and we compute the resulting combinatorial 
quotient 

(.pi = ( 4 ) 

Of course not all elements of T come from geometric configurations. So the task 
is to sort out those which cannot be obtained geometrically and then to find the 
minimum of g(T)’s for those which can. 


3 Criteria for the nonexistence of a geometric configuration 

In m we introduced a number of criteria to deal with this problem. Here we begin 
with a useful modification of what was called a two pencils criterion. We keep this 
name and hope that this will not lead to any confusion. 

Lemma 3.1 (Two pencils criterion). Let L = {L±,..., L,i} be a configuration of 
lines in the projective plane P 2 (K), with the singular set {Pi,... ,P S }, with s ^ 2 . 
Let m\ ,..., m s be the multiplicities of points Pi,... ,P S respectively. Without loss 
of generality we can assume that 

mi ^ m 2 ^ ^ m s . 
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Then either 


m\mi + 2 ^ s 


(5) 


or, if © does not hold, 


(mi - 1 )(m 2 — 1 ) + a ^ s, 


where a is equal to the minimal number of singular points lying on a line passing 
through P\ and P 2 . The number a can be easily computed combinatorially. 

Proof. Assume that points Pi and P 2 do not lie on a configuration lines, then the 
lines from two pencils (lines through Pi, resp. P 2 ) meet in m\m 2 points. Together 
with Pi and P 2 we get ©. 

If © does not hold, points Pi and P 2 lie on a configuration line. Lines from 
these two pencils (apart from the common line) meet in (mi — l)(m 2 — 1) points. 
Now we add the number of points on the common line, which is at least a. □ 

The following Example illustrates how the two pencil criterion is applied 

Example 3.2. The following data: d = 10, t% = 7 and £4 = 4 is a solution of 
©• Then mi = m 2 = 4 and s = 11. Since 4-4 + 2 > 11, we pass to the second 
inequality. Now a = 3, since the line through Pi and P 2 meets with six other lines 
at these two points, hence there must be another point on this line (and one point 
of multiplicity 4 suffices). The inequality 3-3 + 3 > 11 shows that there is no 
geometrical configuration satisfying above data. 

The next idea is to doubly count the incidences. First we need to introduce some 
notation. To a configuration line L we attach its type vector 


u(L) = (z+(L), 1 / 3 (X), ■ ■ •, v d (L)), 


where i^(L) denotes the number of points of multiplicity k on L. For example the 
line L in Figure [T] has type v(L) = (1,2,0,0,0). 



Figure 1 


Now, let n u (L) be the number of lines in L with the type vector u. Then we 
have 



( 6 ) 


X) n u (L) ■ v k = k ■ tk, for k = 2,... ,d. 

u=(v 2 ,...,v d ) 


The first equation simply counts all lines in a configuration. The others count all 
“incidences” — a line passing through a point of given multiplicity k count as one 
incidence. 
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Let T = {i 2 5 ■ • • ,td} be a set of integers satisfying (J2J) for a fixed d. To these 
numbers there is the associated system of equations ([ 6 |). In this system the symbols 
n v (L) are unknown. Which type vectors v(L) can appear in the given configuration 
can be easily determined in advance. Their number is quite restricted. If the system 
© has no-negative solutions, then it follows that the set T cannot be realized 
geometrically. 

However the set of equalities ([6]) is not always sufficient for our purposes. Let 
us assume that in a configuration there is a unique point of multiplicity m, that is 
t m = 1. Then all lines passing through this point (that is all lines L with v(L) = 
(z/ 2 , • • •, I'm- 1 ) 1) t'm+ii • • • i v d)) belong to the same pencil. Now for k ^ m we count 
all points of multiplicity k on these lines, which obviously must be at most tf 


£ 


n, 


,(£) • z/% ^ tk, for k = 2,..., d, k ^ m. 


(7) 


v={v2,—,v d ),Vm=l 


So the new and powerful criterion for nonexistence works as follows: write down 
all equations ([6]) with the set of inequalities flU) for all t m = 1, then try to solve 
this system of linear equations and inequalities in non-negative integers. This is a 
problem, well-known as integer programming, and there are many algorithms and 
software to deal with it. 


Example 3.3. The following data: d = 14, t% = 7, = 10 and t$ = 1 is a solution of 

(J2|) . If it corresponds to a geometrical configuration then there are exactly four type 
vectors v , for which n v { L ) may be non-zero (a fixed line must meet with 13 other 
in singular points using only multiplicities appearing in the configuration, hence 
we can easily write down all possibilities). These are (0, 5,1,0,...), (0, 2, 3,0,...), 

(0,3,1,1,...), (0, 0,3,1,... ). Assume that we have a (resp. b, c, d) lines in L with 
resp. types. We have the following system of equalities: 

a + b + c + d = 14, 

5a + 2b + 3c = 21, 
a + 36 + c + 3d = 40, 
c + d = 5. 

There are two nonnegative integer solutions, namely (a, 6, c, d) € {(0,9,1,4), (1,8, 0, 5)}. 
Observe however that = 1 allows us to use two additional inequalities 

f 3c < 7, 

|c + 3d ^ 10. 

The last inequality gives a contradiction, hence the initial data in this example does 
not come from any geometrical configuration. 


3.1 A Singular script 

We wrote a Singular script, which, given a number of lines d and a value h := q(T ) 
for a geometric configuration, works as follows: 

• it enumerates first all possible arrays of integers T = (< 2 , ■.. ,td) satisfying fl2J), 

• for each such an array it computes the quotient q{T ) as in dU), 



7 


• for those quotients, which satisfy q(T) < h it checks whether the two pencils 
criterion works, 

• if this is not the case, then the script produces an input for linear programming 
problem given by d6|) and fl7|), then it uses a glpsol software to solve it, 

• the results are reported; if for all T with q(T) < h one of the two above criteria 
verifies the non-existence of a geometric configuration with data T, then h is 
a lower bound for H(d). Otherwise the test fails and we do not know H(d) 

The script is revoked by the command 

check(number_of_lines, tested_bound, "output_file"); 


For example check(10,-29/12,”result”) checks the validity of the number 77(10) pro¬ 
vided in Theorem PI 


4 Proofs of the lower and upper bounds 

In this section we prove Theorem 11.61 and Theorem P We begin with the lower 
bound. 


Proof of Theorem 17.61 For d ^ 6 and s)lwe consider the following function 


f{d,s) 


d 1 1 f~ 4 d 2 - 4 d 

~s~2~2\ 1+ s 


For a fixed held K and positive integer d, let £ be a configuration of d lines with 
altogether s singular points. Then we have the following 

Claim 

H(K,£)>f(d,s). (8) 

Taking this for granted, Theorem 11.61 follows easily. Indeed, first of all the right 
hand side in © does not depend on IK, so that 


H(d) ^ min f(d, s ). 

S>1 


If £ is a pencil, i.e. s = 1, then H( IK, T) = f(d, 1) = 0. Otherwise by the celebrated 
de Bruijn-Erdos Theorem [5] it must be s ^ d. Elementary calculus shows that for 
a fixed d the function f(d, s ) is strictly increasing for s ^ d. Hence finally 

H(d) ^ min f(d, s ) = -\y/ 4d - 3 + 

s^d 2 2 


The extra assertion of Theorem P will be proved at the end of this section. 
Now we turn back to the Claim. 

Using © we have 


H(K, L) 


d _ J2l=i m c{Pk) 
s s 


so it suffices to show that 


M:= i i n 

2 2 


4d 2 - 4 d 


s 


s 


(9) 
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The idea now is to apply Jensen’s inequality (fTUj) to (|2|). 

Recall that for a convex function <p(x) and non-negative numbers Ai,..., X s such 
that 1 Aj = 1 there is 

y a i<p{xi) > y a ^ j. (io) 

»= 1 V i=l / 

The function <£>(x) = x(x — 1) satisfies the assumptions. Hence, from Q we obtain 
with Ai = ... = A s = A 

\d{d- 1) = \Ylk=i m d p k)(rn c (P k ) - l) > 

> (J EjLi EjUl m^Pk) - 1) (11) 

= M (M — 1). 

It is elementary to check that this implies © and we are done. □ 

Now we prove the upper bound. 

Proof of Theorem This bound is obtained in a rather naive way. Let r = r(d). 
We consider the projective plane P 2 (F r ) as embedded in the projective plane defined 
over the algebraic closure F r . Then L \ is the configuration of all d\ — r 2 + r +1 lines 
coming from P 2 (F r ). Then we take d 2 = d — d\ general lines in P 2 (F r ). These lines 
form another configuration Since they are general, they intersect pairwise in 
(f l f ) distinct points and they intersect the lines in L i in cfodi distinct points. Thus 
for L = U L 2 w e have 

f \d 2 {d 2 — 1) + did 2 for k = 2 
tk(L) = < r 2 + r + 1 for k = r + 1 . (12) 

[ 0 otherwise 

The bound follows then computing the Harbourne constants and expressing every¬ 
thing in terms of d and r = r(d). Note that for d = r 2 + r +1, we get H(d) ^ —r{d). 

□ 

We conclude this section by showing the extra claim in Theorem 11.61 
Corollary 4.1. Let d = q 2 + q + 1 with q G Q. Then 

H{d) = -q 

and H(d) is computed by the configuration of all lines in P 2 (F g ). 

Proof. The inequality H{d) ^ —q follows from Theorem 11.71 The lower bound 
H(d) ^ — q follows in turn from Theorem 11.61 Note additionally that the proof of 
Theorem 11.61 shows then that there is the equality in ©. Hence d = s in this case 
and we conclude again by the de Bruijn-Erdos Theorem. Note that whereas the 
configuration consists of all all lines in P 2 (F g ), it might be embedded in some larger 
projective plane. □ 
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5 Results justifying Conjecture 11.21 


We show first that there are infinitely many values of d such that 

H(d) ^ h(d) 


holds. More precisely we have the following result 
Proposition 5.1. Let d be a positive integer such that 

{q - l) 2 + {q - 1) + 1 < d < q 2 + q + 1 


with q = q(d). Then 


H(d) ^ h{d). 


Proof. The idea is to construct a configuration of lines with invariants indicated in 
Conjecture 11.21 To this end let % = q 2 + q + 1 — d. Let Lq be the configuration of all 
lines in P 2 (F g ). For i ^ q + 1, we fix a point P\ € P 2 (F 9 ) and remove exactly i lines 
passing through the point P\ getting the configuration These lines intersect only 
at Pi, so that with every line we decrease the number of (q + l)-fold points by q and 
increase the number of g-fold points by q as well. The multiplicity mi of the point 
Pi is q + 1 — i, whereas for i = q and i = < 7 + 1 , the point is no more a singular point 
of the configuration. It is then elementary to check that 

H(¥ q ,L 1 ) = h(d). (13) 

For q + 1 < i ^ 2q — 2, we remove all q + 1 lines passing through Pi. This results in 
a configuration of q 2 lines with q 2 + q points of multiplicity q. Then we remove the 
remaining i — (q + 1) lines from the pencil of lines passing through a second point 
P 2 . Counting as above, we get m 

Finally for i = 2q — 1, we fix two points Pi, P 2 and remove the line joining them, 
and 2 (q — 1) additional lines: q — 1 from a pencil through Pi and q — 1 from the 
other pencil. This results in a configuration of q 2 — q + 2 lines with 


tq+l 

= 1 


tq 

= 3(?-l) 

(14) 

tq—1 

= (.q — i) 2 


and all other tk = 0. This gives H(d) 

^ h{d) also in this case. 

□ 


Now we are in the position to prove Theorem 11.41 

Proof of Theorem ] 1.4\ For all d between 2 and 31, Proposition 15.11 applies, so that 
H (d) is at most equal to the numbers stated in Table |TJ Turning to the lower bound 
it turns out that our Singular script works in all cases. This ends the proof. □ 

We pass now to d in the range 32 ^ d ^ 43. For these values of d we have 
q = q(d) = 7, so consequently d ^ (q — l) 2 + (q — 1) + 1. Hence the construction 
used in Proposition 15.11 does not apply. It is well known that there is no projective 
plane with 43 points (this would correspond to q = 6), see [3] . 

Of course tj = 43 and all other tk = 0 is a solution to m with d = 43. Our 
program cannot exclude this configuration. It also cannot exclude any configuration 
resulting from this fake P 2 (Fg) configuration by removing lines. What it can is to 
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exclude any lower values of Harbourne constants. So that we can conclude that for 
32 ^ d ^ 42 it is 

H(d) > h(d) 

and 

(43) > -6. 

It would be very interesting to modify our approach in a way opening access to 
configurations coming from fake projective planes. We hope to come back to this 
question in the next future. 

Acknowledgement. This paper was to large extend written while the last author 
visited the University of Freiburg. It is a pleasure to thank Stefan Kebekus for 
hospitality. We would like also to thank Tomasz Szemberg for helpful conversations. 


6 Appendix 


ring R=0,x,dp; 

proc writelistH(list 1) { 

string s="H-const: "+string(l[1])+" conf: " 
for (int k=2;k<=size(l);k++) { 
s=s+string(l[k])+" 

} 

return(s); 

} 

proc writelist(list 1) { 
string s=""; 

for (int k=2;k<=size(l);k++) { 
s=s+string(l [k]) + " 

} 

return(s); 

} 

proc computeH(int n, list m) { 
number h=n*n; 

for (int i=2;i<=size(m);i++) { 
h=h-i*i*m [i]; 

> 

int s=0; 

for (i=2;i<=size(m);i++) { 
s=s+m[i] ; 

> 

h=h/s; 
return(h); 

} 

proc twopencil(int n, list m) { 
int s=0; 

for (int i=2;i<=size(m);i++) { 
s=s+m[i] ; 

> 

string result=""; 
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if (s<n) {returnC few points works");} 
int r=2; 
int p=size(m); 
list tp; 
while (r>0) { 
if (m[p]>0) { 
tp[r] =p; 
m [p] =m [p] -1; 
r=r-l; 

} else { 

p=p-i; 

> 

> 

if ((tp [1]-1) * (tp [2]-l)+2>s) -(result=result+ " TP works";} else 

{ 

int dp=2; 

int zp=n-(tp [1]-1)-(tp [2]-1)-1; 
while (zp>0) { 
if (m[p]>0) { 
zp=zp-(p-l); 
m[p]=m [p] -1; 
dp=dp+l; 

} 

else 

1 

p=p-i; 

} 

} 

if ((tp[1]-1)*(tp[2]-l)+dp>s) {result=result+" TP(p) works";} 

} 

return(result); 

} 

proc writeconf(list 1) { 
string s=" line conf: "; 
for (int k=2;k<=size(l);k++) { 
s=s+string(l[k])+" "; 

} 

return(s); 

} 

proc eqcritbyglp(int n, list m) { 

"checking "+writelist(m); 

" conf for n="+string(n)+"..."; 
list rm,nm,sm; 

for (int i=l; i<=size (m); i++) -(rm[i]=0;} 
int p=0; 
int s; 
list cp; 

while (p<=size(m)) { 
rm[2]=rm [2]+1; 

p=2; 
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while ((p<=size(m))&&(rm [p] >m [p])) { 
rm[p]=0; 
p=p+l; 

if (p<=size(m)) {rm[p]=rm[p]+1;} else {break;} 

} 

s=0; 

for (i=2;i<=size(m);i++) { 
s=s+rm[i]*(i-1); 

} 

if (s==n-l) { 

cp[size(cp)+l]=rm; 

} 

} 

int vrb=size(cp); 

if (vrb==0) {returnC C0NF(0) works");} 
list eqs; 
list eqq; 

for (i=l;i<=size(cp);i++) { 
eqq[i]=1; 

} 

eqq[vrb+1] =n; 
eqs[1]=eqq; 

for (p=2;p<=size(m);p++) { 
if (m[p]>0) { 

for (i=l;i<=size(cp);i++) { 
rm=cp [i]; 
eqq[i]=rm[p] ; 

} 

eqq[vrb+1]=p*m[p]; 
eqs[size(eqs)+l]=eqq; 

} 

} 

int j; 

string name=":w test"; 

write(name,"minimize value: al"); 

name=":a test"; 

write(name,"subject to"); 

string wr; 

for (i=l;i<=size(eqs);i++) { 
wr="e"+string(i)+": "; 
eqq=eqs [i]; 

for (j=1;j<=vrb;j++) { 
if (j>l) {wr=wr+" + ";} 
wr=wr+string(eqq[j])+" a"+string(j); 

} 

wr=wr+" = "+string(eqq[vrb+1]); 
write(name,wr); 

} 

int mm,k,o; 

for (mm=2;mm<=size(m);mm++) { 
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if (m [mm] ==1) { 

for (i=2;i<=size(m);i++) { 
if ((i!=mm)&&(m[i]>0)) { 

wr ="b"+string(mm)+"k"+string(i)+": "; 
o=0; 

for (k=l;k<=size(cp);k++) { 
rm=cp [k]; 
if (rm[mm]==l) { 

if (o==0) { 0 = 1 ;}- else {wr=wr+" + 
wr=wr+string(rm[i]) + " a"+string(k); 

} 

> 

wr=wr+" <= "+string(m[i]); 
if (o==l) {write(name,wr);} 

} 

> 

> 

> 

write(name,"integer"); 
for (j=1;j<=vrb;j ++) { 

write(name," a"+string(j)); 

} 

write(name,"end"); 

int dummy=system("sh","glpsol —lp test -o solution"); 
link solfile=":r solution"; 
string sol=read(solfile); 

if (find(sol,"UNDEFINED",1)>10) {return(" SOLVER works");} 
returnC"); 

} 

proc throw(int n, list m) { 
string ii=writelistH(m); 
ii=ii+twopencil(n,m); 

if (find(ii,"works")==0) {ii=ii+eqcritbyglp(n,m);} 
return(ii); 

} 

proc scheck(int n, number bnd, string infofile, list m) { 
int ntp,ntpp,nsolver; 
int ok=l; 

infofile=":a "+infofile; 

write(infofile,"Input data: "+string(n)+" lines, bound for H-constant 
"+string(bnd)+"."); 

write(infofile,"Configurations to exclude:"); 
list v; 
string info; 
int na; 

for (int i=2;i<=n;i++) { 
v[i]=i*(i-l) div 2; 

} 

int sum=n*(n-l) div 2; 
list b; 
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for (i=2;i<=n;i++) { 
b[i]=sum div v[i] ; 

> 

int p; 
int mm; 
while (p<n) { 
m[3]=m[3]+1; 

P=3; 

while ((m [p]>b [p])||((p>n+l-mm)&&(p<mm))) { 
m[p]=0; 
p=p+l; 

m[p]=m [p]+l; 
if (p>mm) {mm=p;} 

> 

if (p==n) {break;} 
m[2]=sum; 

for (i=3;i<=n-l;i++) { 
m [2] =m [2] -v [i] *m [i] ; 

} 

if ((m[2]>=0)&&(m[n]==0)) { 
na=na+l; 

if ((na mod 10000)==0) {string(na) + " already checked..."jwritelist(m);} 
m[l]=computeH(n,m); 
if (m[l]<bnd) { 
info=throw(n,m); 
inf o; 

if (find(info,"works")==0) {ok=0;} 
write(infofile,info); 
if (find(info,"TP ")>0) {ntp++;} 
if (find(info,"TP(p)")>0) {ntpp++;} 
if (find(info,"S0LVER")>0) {nsolver++;} 

} 

} else { 
p=3; 

while (m[p]==0) {p=p+l;} 
m [p] =b [p] ; 

} 

} 

if (ok==l) { 

"All configurations have been excluded."; 

write(infofile,"All configurations have been excluded."); 

> 

info="TP used "+string(ntp)+" times, TP(p) used "+string(ntpp)+" times, 
SOLVER used "+string(nsolver)+" times."; 
inf o; 

write(infofile,info); 

> 

proc check(int n, number bnd, string infofile) { 
list m; 

for (int i=l;i<=n;i++) { 
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m[i] =0; 

> 

scheck(n,bnd,infofile,m); 

} 

proc contcheck(int n, number bnd, string infofile, list sm) { 
list m; 

for (int i=l;i<=n;i++) { 
m[i] =0; 

> 

for (i=l;i<=size(sm);i++) { 
m[i]=sm [i] ; 

} 

scheck(n,bnd,infofile,m); 
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